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SUMMARY 


The longitudinal linearized equations of motion in wind shear have been 
derived for the NASA Terminal Configured Vehicle, a modified Boeing 737 air- 
plane. In addition to the apparent acceleration terms due to wind shear, the 
equations include altitude-dependent stability derivatives. A linear analysis 
of these equations indicates a first-order divergence type of instability due 
to wind shear in which head wind decreases with altitude. Furthermore, this 
instability cannot be stabilized by attitude control alone. However, attitude 
control used in combination with an additional feedback loop which consists 
of the energy height rate feedback to the throttle proved to be effective in 
suppressing instability due to wind shear. In the present report the term 
"energy height" denotes the sum of kinetic and potential energies per unit 
weight of the airplane referenced to the surrounding air mass. A brief 
piloted, real-time, nonlinear simulation indicated the desirability of using 
a display based on the rate of change of energy height and of commanded thrust. 


INTRODUCTION 

The effect of wind shear on aircraft trajectories and aircraft control 
has been the subject of numerous studies dealing with aviation safety. These 
studies are usually conducted on complex simulations of specific aircraft with 
an assumed atmospheric wind model which includes a particular variation of the 
wind as a function of Earth-fixed coordinates. In contrast, references 1 and 2 
deal with the effect of wind shear on the longitudinal stability and control 
within the framework of a small-disturbance theory of flight dynamics. The 
analysis presented here is an extension of the method proposed in reference 1, 
in that it includes stability derivatives with respect to height in addition 
to the kinematic effects of wind shear. 

In the mathematical model chosen for the study, the parameters were those 
of the Boeing 737-100 airplane currently used in the NASA Terminal Configured 
Vehicle (TCV) Program, A linear variation of the wind magnitude with height 
was assumed in the analysis. This assumption may limit the applicability of 
the results, since actual flight experience indicates that not only the magni- 
tude, but also the direction, of the wind may change with altitude. Such 
changes are also frequently accompanied by sustained updrafts and downdrafts. 
Although a rigorous treatment of the effect of these types of atmospheric 
motion on longitudinal stability is beyond the scope of the present analysis, 
the control law proposed here for alleviating the effects of linear wind shear 
is shown to be useful in alleviating the effects of updrafts and downdrafts. 


SYMBOLS 

Values are given in both SI and U.S. Customary Units. Measurements were 
made in U.S. Customary Units. 


A coefficient matrix in linear equations of motion 

B control matrix in linear equations of motion 

C L lift coefficient 

C L « nondimensional rotary stability derivative, ?C h /'d(c'a/2V) 

C m pitching-moment coefficient 

Cma nondimensional rotary stability derivative, 9Cn/9(cc/2V) 

c mean aerodynamic chord, 3.41 m (11.2 ft) 

D aerodynamic drag, N (lbf) 



3d 

D v 

= — / kg/sec (slugs/sec) 
oV 

D a 

9d 

= ~ t kg-m/sec 2 (slug-f t/sec 2 ) 
3a 


3d 

°6 e 

= — , N/deg (lbf /deg) 
3o e 


3D 

d 6t 

= 

3(6t) 


F X ' F z total force components in the direction of x- and z-wind axes, 
respectively, N (lbf) 

g gravitational acceleration, 9.81 m/sec 2 (32.2 ft/sec 2 ) 

h altitude above mean sea level, m (ft) 

, V 2 

h e energy height, — + h, m (ft) 

2g 

Im(s) imaginary part of complex variable 

Iyy moment of inertia, kg-m 2 (slug-ft 2 ) 

K h e proportional gain, deg/ (m/sec) (deg/(f t/sec) ) 

Kl integral gain, deg/m (deg/ft) 
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K 6 gain on pitch rate to elevator, deg/deg 

L aerodynamic lift, N (lbf) 


L q 

3l 
= 3q # 

N-sec (lbf-sec) 

w 

3l 
“ 3v' 

kg/sec (slug s/sec ) 

% 

3L 
" 3a' 

kg-m/sec 2 (slug-ft/sec 2 ) 


3l 

~ TT / 

3a 

kg-m/sec (slug-ft/sec) 

L «e 

3l 

= 86? 

N/deg (lbf /deg) 

M 

aerodynamic pitching moment, N-m (lbf-ft) 

Mg 

3M 

= 3q' 

kg-m 2 /sec (slug-ft 2 /sec) 

M v 

3m 
~ 3v' 

kg-m/sec (slug-ft/sec) 

% 

3m 
” 3a' 

kg-m 2 /sec 2 (slug-ft 2 /sec 2 ) 

M«i 

3M 

3a' 

kg-m 2 /sec (slug-ft 2 /sec) 

M S e 

3m 
= ™~e 

kg-m 2 / (sec 2 -deg) (slug-ft 2 /(sec 2 -deg) ) 

m 6t 

3M 

-, m (ft) 


Ui 3 (<$T) 

m mass, kg (slugs) 
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q pitching velocity about body Y-axis, rad/sec 

q w pitching velocity about wind y-axis, rad/sec 

Re(s) real part of complex variable 

S reference area, m 2 (ft 2 ) 

s variable in Laplace transformation, sec - ! 

T total thrust, N (lbf) 

T e trim value of thrust, N (lbf) 

3T 

T V = —r kg/sec (slugs/sec) 

oV 

t time, sec 

u w magnitude of horizontal wind, m/sec (ft/sec) 

du w 

' w -1 

u w = , sec 1 

dZ E 

V total velocity relative to moving air mass, m/sec (knots) 

V e equilibrium velocity relative to moving air mass, m/sec (knots) 

v w magnitude of lateral wind, m/sec (ft/sec) 

w w magnitude of vertical wind, m/sec (ft/sec) 

Z E Earth-fixed Z-coordinate, m (ft) 

a angle between zero lift line and relative wind vector, rad 

ol g trim value of angle of attack, rad 

dr angle between thrust line and zero lift line, rad 

y angle between local horizontal and relative wind, rad 

y e equilibrium or trim value of y, rad 

Ah perturbation in altitude from the initial condition, m (ft) 

AV perturbation in relative wind, m/sec (ft/sec) 

AZ e perturbation in Z E -coordinate, m (ft) 
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Aa 

Ay 

A0 

6 e 

5 e,p 

6 e,STAB 

6T 

6 th 

^th,c 

0 

0 ref 

0 w 

£p' ^sp 
W P» W sp 

A 


perturbation angle of attack, rad 

perturbation in relative flight-path angle, rad 

perturbation in pitch attitude, rad 

elevator deflection measured from trim value, deg 

elevator command from pilot, deg 

trim value of stabilizer position, deg 

thrust perturbation measured from trim value, N (lbf) 

throttle position, deg 

throttle position command, deg 

pitch attitude, deg 

reference value of pitch attitude, deg 

Euler angle of x-wind axis measured from Earth-fixed X-axis, deg 

damping ratio of phugoid and short-period modes, respectively 

undamped natural frequency of phugoid and short-period modes, 
respectively, rad/sec 

dot over a symbol denotes the first derivative taken with respect to 


time. 


MODEL DEFINITION 

In this section, the wind is defined by the same simple linear relation- 
ship which was used in reference 1 : 

u w = u w z E ( 


du. 


where u w = 


dZ 


E 


v„ = 0 


w w = ° 


(2) 

(3) 
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That is, the only nonzero component of the wind lies in the horizontal plane, 
parallel to the Earth-fixed X-axis, and this component varies linearly with the 
Earth-fixed Z-ooordinate Zg. By definition, u w is positive in the case of 
a tail wind. In all cases the wind gradient will be considered a constant. 

Figure 1 illustrates the definition of the axis systems and the wind gradient. 

The airplane is represented by a longitudinal small-disturbance model 
wherein wind axes are used for the lift-force and drag-force equations, and body 
axes for the pitching-moment equation. The small-disturbance model is obtained 
by linearizing about an equilibrium path the following nonlinear equations of 
motion: 


* | 

m(V - Vu w sin 9 W cos 9 W ) = F x ( 4 ) 

m (“Vq w - Vu^ sin 2 9 W ) = F z ( 5 ) 


I yy q = M 


(6) 


The derivation of these equations is presented, with minor differences in nota- 
tion, in reference 1 . 

The equilibrium flight path is taken to be a path along which V, q w , 
and q are equal to zero. Since 


<3w = q - « 


(7) 


the steady-state flight path is seen to be either level flight, steady climb, 
or dive with fixed controls. Although this flight path corresponds to generally 
accepted notions about flight in a steady state, reference 1 shows that the 
definition just given of steady state allows flight paths which do not appear 
rectilinear for an observer fixed to an inertial reference frame. 

Except for the terms due to the wind shear, equations (4), (5), and (6) 
are those given in standard textbooks on airplane stability and control. Con- 
sequently, the reader is referred to reference 3 (pp. 157-163) for details of 
the linearization. The additional acceleration terms which appear on the left- 
hand side of equations (4) and (5) can be easily linearized. Since this report 
deals with two-dimensional motion, the symbol 3 W denoting the flight-path 
angle is replaced by the more conventional symbol y. A Taylor series expan- 
sion to the first order results in the following expressions of the terms which 
contain the wind gradient: 
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Vu^, sin y cos y = V e u^ sin y e cos y e + ^- sin 2y e j A V 

+ (V e u^ cos 2y e ) Ay (8) 

Vu^ sin 2 y = V e u^ sin 2 y e + (u^ sin 2 y e ) AV + (V e u^ sin 2y e ) Ay (9) 

where the higher order terms have been omitted and the following convention has 
been used: 


V = v e + AV 

(10) 

y = y e + Ay 

(11) 


Thus, the presence of wind shear modifies Etkin's linearization by the addition 
of the equilibrium or trim terms 


I 

-V e u„ sin y e cos y e 


and 


-V e u^ sin 2 y e 


to the equations describing the trimmed flight. Also, the linearized drag-force 
and lift-force equations of Etkin (ref. 3, p. 163) must include the following 
terms which result from wind shear: 

Drag-force equation: 


- u w sin 2y e ^ 


AV - (V e u^ cos 2y e ) Act + (V e u^ cos 2y e ) A9 


Lift-force equation: 


(u^ sin 2 y e ) AV - (V e u^ sin 2y e ) Aa + (V e u^ sin 2y e ) 


Ae 
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Stability Derivatives Due to Altitude Change 

These derivatives are normally neglected because most linear analyses 
assume the atmosphere to be uniform. In the present analysis, however, these 
Zg-derivatives naturally arise since the wind speed varies with the 
Zg-coordinate. However, the Zg-der i vatives will be included in the analysis 
only for the case of horizontal, wings-level, trimmed flight; i.e., the analysis 
will consider only two cases: 


1 . Y e = 0 and Z E -deri vatives are nonzero 


The 


Z E 


2. y e ^ 0 and Zg-derivatives are assumed to be zero 


8F X 

9F Z 

3m 

derivatives , 

, and 

are evaluated as 

Zg 

Zg 

az E 


3f x 

3f x 

3uw 

3F X 

1 

3z e 

9u w 

3z e 

~ 3v~ 

U W 

3f z 

3F Z 

3u w 

3f z 

i 

3z e 

3u w 

3z e 

~ 3V~ 

“w 

3M 

3m 

3u w 

3m 

1 

3Zg 

3u w 

3Zg 

” 3v" 

u w 


( 12 ) 


( 13 ) 


( 14 ) 


This is, for horizontal flight, the Zg-derivatives can be obtained very simply 
from the conventional speed-derivatives if the wind gradient u^ is known. For 
diving or climbing flight, the evaluation of these derivatives would be much 
more complicated because of their dependence on altitude, flight-path angle, 
atmospheric density, etc. 

Thus, the longitudinal equations of motion of an airplane flying in a ver- 
tical wind gradient can be composed fran the conventional small-disturbance 
equations by adding certain kinematic terms and the Zg-derivatives. For the 
sake of completeness, both the conventional and wind shear terms are given in 
their literal form in the appendix. 


Linearized Equations of Motion of the TCV Airplane 

The specific airplane selected for the analysis is the NASA Terminal 
Configured Vehicle (TCV), which is a modified Boeing 737. The linear equations 
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of motion were evaluated in level flight and in a 3° descent at the flight 
conditions described in table I (vehicle parameters are also described). 


TABLE I.- FLIGHT CONDITIONS AND VEHICLE PARAMETERS 


Altitude, m (ft) 

Equivalent airspeed, m/sec (knots) 

Weight, kN (lbf) 

Iyy, kg-m 2 (slug-ft 2 ) 

Gear 

Flaps, deg 

Reference area, m 2 (ft 2 ) 

Mean aerodynamic chord c, m (ft) 
Center of gravity, percent c . . 


152.4 (500) 

64.77 (125) 

378 (85 000) 

1.08 x 10 6 (7.99 x 10 5 ) 

Down 

40 

91.04 (980) 

3.41 (11.2) 

18 


The coefficient matrices of the linearized equations of motion are listed 
in the appendix. The matrices were computed for the no-shear case by using the 
computer program described in reference 4. Wind shear effects were included by 
adding both the kinematic terms and the Z^-derivatives to the computed matrices. 
In table II trim values of the more important vehicle parameters are listed. 


TABLE II.- TRIM VALUES OF THE TCV AIRPLANE PARAMETERS 


Parameter 

0 

o 

II 

<U 

>- 

n> 

ii 

i 

u> 

0 

V e , m/sec (knots) 

64.77 (125) 

64.77 (125) 

a e / de 9 

2.28 

2.38 

T e , kN (lbf) 

57.6 (13 174) 

38.1 (8712) 

5 e,STAB/ de 9 

8.37 

8.34 

6 e , deg 

2.76 

2.54 


EFFECT OF WIND SHEAR ON THE BASIC AIPRLANE 

The longitudinal motion of the basic airplane in the absence of wind shear 
is characterized by an adequately damped short-period mode (w S p =1.35 rad/sec; 
C S p = 0.46) and a lightly damped phugoid mode (Wp = 0.166 rad/sec; £p = 0.08). 
As shown in reference 1 , the principal effect of the type of wind shear assumed 
here is reflected in the phugoid motion of conventional aircraft. This was 
certainly the case for the TCV airplane. Varying the wind gradient u^ between 
-0.1 sec - ! and 0.1 sec -1 had an insignificant effect on the short-period mode. 
The effect of varying Uw on the phugoid characteristics is depicted in the 
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root contours shown in figure 2 for y e = 0° and y e = -3°. There is a sig- 
nificant difference between the two contours, especially for positive values 
of u^. If u^ exceeds approximately 0.07, the oscillatory phugoid root splits 
into a subsidence and a divergence or "tuck” mode. An additional pole, due to 
the coupled Z E -equation in the y e = 0° case, remains close to the origin at 
all values of u^; hence, this pole is not shown. 

The difference between the two contours (shown in fig. 2) may be attrib- 
uted to either of two factors: the different values of y e or the presence 

or absence of the Z E -der ivatives . To resolve this question, the root contour 
was computed for the y e = 0 case, but with all Z E -derivatives set equal to 
zero. The result is shown in figure 3 along with the original y e = 0 contour 
reproduced from figure 2. The influence of the Z E -der ivatives on the equations 
of motion is much greater than the influence of y e . In the computation of the 
poles for the y e = 0 case in the presence of wind shear, the following obser- 
vation (pointed out by Windsor L. Sherman, of the Langley Research Center) can 
be made: The shift in the location of the poles in the complex plane is the 

same whether the shift is due to the Z E -der ivatives or to the additional accel- 
eration terms in equations (8) and (9) . 

For purposes of comparison, time histories were computed for the Y e = 0 
cases with three values of the wind gradient u w . The values selected were 
u w = 0, i.e., no shear, and u w = ±1 sec” 1 2 3 . The latter values correspond to 
very large shear, and it is unlikely that wind shears of this magnitude occur 
in an atmospheric layer thick enough for the airplane to achieve equilibrium. 
These extremes were chosen to demonstrate clearly the effect of wind shear on 
the motion. Figure 4 illustrates the controls-f ixed response of the airplane 
to an initial step tail wind gust of 3 m/sec. In comparison with the no-shear 
(u^ = 0) case, and the case in which head wind increases with decreasing alti- 
tude (u w = -1), the traces showing the motion in a wind shear when u w = 1 
reveal the following information about the latter type of shear. Up to about 
14 sec into the time histories, shear of this type is accompanied by both an 
airspeed and an altitude loss caused by the initial disturbance. The u w = 1 
case clearly shows the divergence or "tuck" mode due to wind shear. It should 
be pointed out that the time histories were computed with the controls fixed; 
normally, the pilot would not allow the large speed or altitude excursions to 
occur . 


Attitude Control in Wind Shear 

The simplest form of attitude control is the feedback of pitch attitude 0 
to the elevator S e through a pure gain Kg, as shown schematically in fig- 
ure 5. The root locii presented in figure 6 depict the closed-loop character- 
istics of this type of control under the same three types of wind shear con- 
ditions for which the time histories were computed in the previous section: 

1 . Head wind decreasing with altitude gain, u^ = -1 

2. No wind shear, u^ = 0 

3. Head wind increasing with altitude gain, u w = 1 
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All three root locii show qualitatively that the 0 to 6 e feedback does 
not add additional damping to the system; it merely redistributes the total 
available damping between the short-period and the phugoid modes, as reflected 
in the fact that the sum of all the poles does not vary significantly as the 
gain is changed. Figure 6(c) shows the existence of a nonrainimum phase zero in 
wind shear where = 1 . This means that an instability will exist for any 
value of the gain Kq so that neither the pilot nor an attitude feedback con- 
trol system could stabilize the first-order divergent mode by pitch attitude 
closure alone. 


Energy Considerations 

The energy balance of the airplane is affected by the variation of hori- 
zontal wind. This can be shown by multiplying the drag equation by V and 
noting that h = V sin y. The resulting energy relationship 


i(i 

dt\2g 


V 2 




cos a 



u w sin y cos y 


shows that the rate of energy change per unit weight computed relative to the 
surrounding air mass, sometimes denoted by h e , or energy height, is equal to 

V 

the specific excess power — (T cos a - D) plus a term due to the wind shear. 

m 9 y2 ^ 

Note that the wind shear term — u^ sin y cos y can take either a positive 

9 

or a negative sign. During a landing approach, when y < 0, the shear in which 
u^ > 0 results in an excessive decrease of energy height which is due to wind 
shear only. The loss of energy is further aggravated by the ineffectiveness of 
attitude control in stabilizing the first-order divergence (as shown in the pre- 
ceding section) . 

Examination of the foregoing energy equation suggests a possible method 
to cancel the wind shear term as reported by Joppa in reference 5. If specific 
excess power is assumed constant, the presence of wind shear would be indi- 
cated by a change in the energy height rate. This also means that it should 
be possible to vary the thrust term in such a way as to keep the energy rate 
constant regardless of the sign and magnitude of the wind shear term. At con- 
stant airspeed, of course, the energy height rate is just the rate of climb/ 
descent. If approach and landing take place while the airplane is flying on 
the back side of the drag polar, pilots generally control rate descent by 
adjusting the throttle setting. Thus, the closure of energy height rate to 
throttle appears to be a natural choice from the standpoint of both the energy 
balance of the airplane and manual control techniques during landing. Refer- 
ence 6 describes a practical, inexpensive method of airborne energy height 
measurement which is the subject of research studies at the Langley Research 
Center . 
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Closed-loop Control of Energy Height Rate 


In this section, it is assumed that the pilot controls pitch attitude by 
the elevator; thus, the pitch attitude to elevator control loop remained closed 
in the manner previously described with Kq = 0.32. While this choice was some- 
what arbitrary, it resulted in a satisfactory short-period damping ratio of 
approximately 0.4 for all three wind gradient cases: 


u w = "I 
u w = 0 
u w = 1 


For the first two cases, this choice of Kq also increased the phugoid damping, 
but the u^ = 1 case still had the unstable divergent mode. Stabilization of 
the latter was attempted by closing the energy height rate throttle loop through 
the following transfer function: 


6 th,c 

0.5 

6 th 

N 

2651 

6T 


deg 

s + 0. 5 


deg 




The preceding linear gain-lag combination represents an adequate approximation 
of the TCV airplane power plant at the flight condition specified in table II. 
Figure 7 shows the closed— loop configuration examined in the subsequent para- 
graphs. The configuration includes a forward-loop integration which was 
intended to improve the low-frequency characteristics of the h e loop. Fig- 
ure 8 shows the root locii for the three values of u*J without the forward- 
loop integration. Energy height rate feedback is seen to stabilize the unstable 
mode due to wind shear while leaving the short-period characteristics virtually 
unchanged. The only adverse effect of h e feedback is on the phugoid mode at 
u„ = -1. As the gain K£ e approaches infinity, the phugoid pole closes on a 


zero in the right-hand half-plane. Figure 9 shows the same three root locii 

with a forward-loop integration and K x equal to - Kfc . For moderate values 

5 e 

of the gain, say 0.4 < KjJ < 0.5, both the augmented short-period and phugoid 


modes possess acceptable characteristics with or without forward-loop integra- 
tion. Time histories of several motion variables as well as those of the 
elevator and throttle are presented in figures 10 and 11 for the closed-loop 
case Kj = 0.1 and K£ e =0.5. 
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Piloted Simulation of Energy Height Rate Control 


In the preceding linear analysis, the effectiveness of h e feedback was 
shown in stabilizing the divergent mode in the wind shear in which u^J = 1 . 

In order to show that this feedback is also effective in controlling the flight 
path during nonlevel flight and in a more realistic wind shear environment, a 
brief piloted simulation was conducted using a large-scale, nonlinear simula- 
tion of the TCV airplane and the aft flight deck. For a detailed description 
of this simulation, see reference 7. 

In the simulation, the pilot's task was to make an instrument approach on 
a 3° glide slope using either an experimental attitude control wheel steering 
mode or a fully coupled automatic landing mode to apply elevator commands. 
Simultaneously, he controlled engine thrust by making throttle commands either 
through a manual, proportional system or through an engine height rate control 
system in which power lever position was proportional to energy height rate. 

A simplified schematic of the control modes is shown in figure 12. The 
details of the elevator control modes are not shown because the objective of 
the simulation was the comparison of manual throttle control with the energy 
height rate system from the standpoint of pilot work load during wind shear 
conditions. The simulated wind profile (shown in fig. 13) included not only 
a variation of the horizontal wind with altitude, but also two separate regions 
of strong downdrafts. 

A typical approach using the autoland control mode for elevator control with 
manual control of the throttles is shown in figure 14. The task consisted of a 
straight-in approach in approximately level flight until the glide slope was 
intercepted at an altitude of approximately 300 m (1000 ft). Thereafter, the 
pilot attempted to maintain airspeed at 120 knots by manipulating the throttles. 
The magnitude of h e was not displayed to the pilot. The time histories show 
the coupling existing between the power setting and the energy height rate. The 
minimum value of the energy height rate occurred at about 42 sec, approximately 
150 m (490 ft) altitude, as a result of reducing the throttle setting about 
10 sec earlier. Noting the loss of airspeed, the pilot then rapidly advanced 
the throttles to almost full power. The autoland mode almost simultaneously 
commanded a nose-down elevator of about 4°. Despite the nose-down command the 
airspeed continued to fall, and the rate of descent increased, because of the 
decreasing headwind and the increasing downdraft for about 20 sec. The energy 
height rate trace also reflects the combined effect of the horizontal and 
vertical wind. Subsequently, as the downdraft abated and the throttle setting 
remained very high, the energy height rate reached its maximum, and airspeed 
increased significantly. Eventually a relatively smooth touchdown with a 
vertical velocity of 0.6 m/sec (2 ft/sec) was made by the autoland system. 

These results should be compared with figure 15 which depicts time histories 
obtained under the same conditions; the throttles, however, were driven com- 
pletely by the energy height rate system except for manually repositioning the 
throttle levers at the time of glide slope interception and at flare initia- 
tion. (It is significant that during the approach the pilot did not manipulate 
the throttles manually except initially and during flare.) As expected, the 
energy height trace is much smoother than during the previous run, and the 
value of ti e gradually increased from the minimum as the downdraft abated and 
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the horizontal wind became constant. The increase in h e resulted in a 
decreased power setting just before initiating flare. At the expense of the 
airspeed momentarily falling below 105 knots, the large airspeed increase which 
had been experienced during manual throttle operation at this point on the 
approach trajectory was avoided. The approach resulted in a touchdown similar 
to that of the previous case. Detailed comparison of the two approaches indi- 
cates that the use of an energy height rate system in the throttle loop results 
in an approach very similar to and in sane respects slightly better than that 
obtained with manual throttle. The significant point, however, is that pilot 
work load was reduced with the energy height rate system. 

Figure 16 shows an approach flown with attitude-control-wheel steering 
mode in the elevator loop and energy height rate system in the throttle loop. 
This combination resulted in a significantly higher pilot work load than what 
was experienced with the autoland and energy height rate combination. Also, in 
terms of glide slope error and vertical velocity at touchdown, this approach 
trajectory was inferior to those of the previous two cases. The pilot stated 
that the effects of wind shear and downdraft made the pitch axis control diffi- 
cult to recognize as an attitude control system. In the opinion of the pilot, 
manual throttle operation in combination with experimental control wheel steer- 
ing mode would have been completely unacceptable. 

However, despite its beneficial effect in reducing the pilot work load in 
wind shear, the pilot expressed dislike of the throttle loop augmentation by 
the h e signal when used either with the autoland or the experimental attitude 
control modes. The reason stated was that manual throttle commands would be 
canceled or overridden by the h e feedback, resulting in the objectionable 
impression that the engines were running out of control. Excessive throttle 
activity would also preclude the operational use of an automatic ti e to 6T 
loop. 


Consequently, it is postulated that use of the h e signal for energy 
management during operation in wind shear or up/down drafts could best be 
accomplished by displaying the output of the energy height rate device for use 
in manually positioning the throttles. The pilot would then be free to follow 
or disregard its output. A conceptual display for this purpose is depicted 
schematically in figure 17. Typical use of this display is as follows: Prior 

to activating the display, the pilot would null the indicator while flying on 
the desired flight path at the appropriate power setting. After pulling out 
the RESET knob, the pilot keeps the indicator at the null position by appro- 
priately positioning the throttle. As used in this mode, the display serves 
to alert the pilot to the presence of wind shear but leaves him the responsi- 
bility of deciding whether to go around or to continue the approach using the 
output of the display as a guide to modulating the throttle position. Further 
simulation studies are needed to determine both the display format and its 
optimum use during the landing approach. 


CONCLUDING REMARKS 

The longitudinal linearized equations of motion have been derived for the 
Terminal Configured Vehicle (TCV) Boeing 737 airplane and contain not only the 
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apparent acceleration terms due to wind shear, but also some altitude-dependent 
stability derivatives. The inclusion of the latter derivatives amplifies the 
effects of wind shear on the longitudinal modal characteristics. These wind 
shear effects, which have been derived in the general case previously, are also 
true for the TCV airplane. 

In the case of the airplane flying with fixed controls, the type of shear 
in which head wind is decreasing with decreasing altitude results in a first- 
order divergence. Wind shear in which the head wind increases with decreasing 
altitude results in decreased damping ratio and increased frequency of the 
phugoid motion. 

The short-period characteristics remain virtually unaffected in both kinds 
of wind shear . 

The phugoid mode instability due to wind shear cannot be stabilized by 
attitude feedback to the elevator for the TCV airplane; however, adding an 
energy height rate loop to the throttle resulted in a stable phugoid for the 
types of wind shear considered in the analysis. 

A brief piloted simulation was conducted subsequent to the linear analysis, 
largely to verify the results of the latter in a large-scale, nonlinear simu- 
lation of the TCV airplane. The results are of a preliminary nature, but they 
do substantiate the findings of the linear analysis. The piloted simulation 
also indicated that the energy height rate to throttle closure should be made 
through the pilot utilizing a display. Further simulation effort is needed to 
arrive at the most suitable form of this display. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 4, 1980 
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APPENDIX 


LITERAL AND NUMERICAL FORM OF THE COEFFICIENT MATRICES 


The longitudinal small- disturbance equations of motion can be written as 


x = Ax + Bu 

where the variable x is the five-dimensional state vector: 
x = (Av, Aot, q, A0, Az e ) 

and the variable u is the two-dimensional control vector: 
u = (ST, S e ) 


The literal forms of the element of the matrices A and B are listed below, giving 
both the conventional terms (ref. 3) and the terms which must be added if wind 
shear is included in the analysis. The wind shear contribution included the 
newly derived altitude- dependent derivatives (e.g., A] 5, A25, A35) : 


Altitude-dependent 

derivatives 

Conventional term 

A 11 

1 

-{T v cos a T - Dy) 
ra 

A 1 2 

1 

g cos y e (T e sin a T + D a ) 

ra 

A 1 3 

0 

A 1 4 

-9 cos Y e 

A 1 5 

0 

A 21 

Ly + Ty Sin 0^ 
mV e + 


* Mind shear contribution 

\ sin 2 Ye 

-u> e cos 2y e 
0 

u^V e cos 2y e 
1 

(T v cos a T - Dyji^ (zero if y e / °) 

ra 

mu^ sin 2 y e 
mV e + ^ 


a 22 


+ T e cos cup - mg sin Ye 
mV e + 


mu^V e sin 2y e 
mV e + 
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Altitude-dependent 

derivatives 

A 23 

a 24 

a 25 

a 31 

a 32 

a 33 

a 34 

a 35 

a 41 

A42 

a 43 

a 44 

a 45 

a 51 

a 52 

a 53 

a 54 

a 55 


Conventional term 

Lg - mV e 
mV e + Li 

mg sin y e 

mV e + 


Wind shear contribution 

0 

mu> e sin 2y e 
mV e + 1^ 


Lv + Ty sin a T 

(zero if Y e + 0) 

mV e + la 

Mimu^ sin 2 y e 

Iyy(mV e + 1^) 

Mimu^V e sin 2 Y e 

Iyy(mV e + Lij 

0 

>£mu> e sin 2y e 
lyy (mV e + 1^) 


0 0 

0 0 

1 0 

0 0 

0 0 

-sin y e 0 

v e cos y e 0 

0 0 

-V e cos Y e o 

0 0 


/ Ly + Ty sin 

“ : — l M v - Mi : 

iyv V mVp + Ln i 


{ zero if y e / 


1 / Ly + T sin a, 

— My - ^ 1 : 

z yy\ raV e + *a 


•j / La + T e cos a T - mg sin y e ^ 

T — " Ma — 


-lyy 


raV e + i£ 


1 / mV e - Lg 

— Mq + Mi r 

T yy\ mV e + la. 


Mimg sin y e 

Iyy(mV e + I£) 


1 
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Altitude-dependent 

derivatives 


Conventional term 


Wind shear contribution 


B n 

1 

“(cos Op - D 6T ) 

0 


1 


B 1 2 

“ “ B 6 e 
m e 

0 


B21 


b 22 


sin OLj* 
mV e + 

L *e 
mV e + 


B 31 


1 / sin d™ \ 

Kt - T~7 

x yy\ mv e + la/ 


b 32 




— M<S 0 - Mi 


1 yy\ e 


mv e + La 


b 41 b 52 


As mentioned in the section "Model Definition," the conventional derivatives 
were computed by the numerical differentiation technique described in refer- 
ence 2. The terms containing u^, the wind gradient, were computed using the 
above formulas. For that computation, the dimensional rotary stability deriv- 
atives and M^ were needed. These quantities were obtained from their 

defining formulas: 


Lev — qS Cr • 

^ 2 ^ 


(where = -8) 


Mr* ~ qS — ~ C™: 
^ 2V e 


(St = _2 ‘ 65) 


For the linear analysis, the linearized expression for the energy height 
rate was needed. This was approximated by 



v e d 

(AV - 

g dt 


AZg) 
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For the sake of completeness the numerical value of the A-matrix is given for 
values of u^ equal to -1 , 0, and 1 : 


u w = -1 


-0.0459 

11 .67 

0.0 

-16.27 

-0.0046 

.0049 

-.7080 

.9990 

.0 

-.0003 

.0024 

-.4500 

-.529 

.0 

.0003 

.0 

.0 

1 .0000 

.0 

.0 

.0 

64.92 

.0 

-64.92 

.0 



u w = 0 



-0.0459 

5.18 

0.0 

-9.81 

0.0 

-.0049 

-.7080 

.9999 

.0 

.0 

-.0034 

-1 .4500 

-.5290 

.0 

.0 

.0 

.0 

1.0000 

.0 

.0 

.0 

64.92 

.0 

-64.92 

.0 




u w = 1 



-0.0459 

-1 .29 

0.0 

-3.32 

0.0046 

.0049 

-.7080 

.9990 

.0 

.0003 

.0023 

-1 .4500 

-.5290 

.0 

.0003 

.0 

.0 

1 .0000 

.0 

.0 

.0 

64.92 

.0 

-64.92 

.0 

The B-matrix, 

whose elements 

do not depend on 

u^, is also given: 



0.0001 0.0 

.0 -.0008 

.0 -.0213 

.0 .0 

.0 .0 
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Figure 13.- Simulated wind profile. 
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Figure 14,- Time history showing approach flown in wind shear with 

autoland and manual throttle. 
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Figure 16.- Time history showing approach flow in wind shear with attitude 
control-wheel steering in elevator loop and energy height rate feedback 
in throttle loop. 
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